Twinlike models for kinks and compactons in flat and warped spacetime 
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This work deals with the presence of twinlike models in scalar field theories. We show how to build 
distinct scalar field theories having the same extended solution, with the same energy density and the 
very same linear stability. Here, however, we start from a given but generalized scalar field theory, 
and we construct the corresponding twin model, which also engenders generalized dynamics. We 
investigate how the twinlike models arise in flat and in curved spacetime. In the curved spacetime, 
we consider a braneword model with the warp factor controlling the spacetime geometry with a 
single extra dimension of infinite extent. In particular, we study linear stability in the flat and 
curved spacetime, and in the case of curved spacetime, in both the gravity and the scalar field 
sectors, for the two braneworld models. 

PACS numbers: 11.27.+d, ll.10.Kk 



I. INTRODUCTION 

Topological structures are of great interest in high en- 
ergy physics [l|, [H and in other areas of nonlinear science 
[3|-l5|]. In high energy physics, they are kinks, vortices, 
monopolcs and other field configurations. They are fi- 
nite energy field configurations which solve the equations 
of motion of the respective models in the corresponding 
spacetime dimensions. Kinks are the simplest structures, 
and they are usually constructed under the presence of 
real scalar fields in (1,1) spacetime dimensions. 

In this work we focus on the presence of kinks in model 
described by real scalar fields, driving attention to the re- 
cent issue concerning the investigation of twinlike mod- 
els. This issue was firstly considered in [(|, and then 
studied in a diversity of contexts in the works @-[Hj]. 
In the papers on twinlike models 

ji-IHI, the key issue 
was to construct and investigate twinlike models, start- 
ing from a standard model and then introducing the twin 
model, which usually appears describing generalized or 
non-standard dynamics. The important point is that it is 
sometimes possible to find a standard and another model, 
with generalized dynamics, both having the same defect 
solution, with the very same energy density. The two 
models are then twinlike models. However, in Ref. Q 
some of us have shown that it is possible to have twinlike 
models with the very same stability behavior. We call 
this the strong condiction: that is, there are twin mod- 
els, if they have the same solution, with the very same 
energy density; but there are models that are twins in the 
strong sense, if they also have the very same behavior, 
concerning linear stability. 

To enlarge the scope of the work, we also investigate 
models of the Randall-Sundrum type |l2l|. in the presence 
of scalar fields, as suggested in Ref. [l3| . The issue was in- 
vestigated in several works flil . Il5j , and here we consider 
the scalar field with generalized dynamics, leading us to 
a mathematical framework which is much more compli- 
cated then it appears in the case of standard dynamics. 



In spite of this, we introduce a complete investigation of 
linear stability, both in the flat and curved spacetime. 

The main issue here is to open another route to deal 
with twinlike models. Indeed, we consider the construc- 
tion of twinlike models, but now we start from a general- 
ized model, instead of using a standard field theory. The 
issue is of current interest, mainly because models with 
generalized dynamics have been used to respond for the 
presence of dark energy and dark matter, and to test pos- 
sible modifications of general relativity. However, this is 
not a simple question, because of the intricacy of the the 
models to be investigated. To ease the investigation, we 
follow Refs. [lB-O- In particular, we introduce a new 
function W = W((f>), from which we obtain simple first- 
order equations that very much help us to study and solve 
the equations of motion. The presence of W{4>) allows for 
supersymmetric extensions, as previouslly investigated in 

mm 

We organize the current work as follows: in Sec. |TT] 
we introduce the procedure, starting from a generalized 
model. In Sec. IIIII we investigate other generalized mod- 
els, to show that the procedure is general, and works not 
only for the model introduced in the previous section. In 
Sec. IIVI we extend the procedure to the braneworld sce- 
nario, in the case of warped geometry with a single extra 
dimension of infinite extent. Finally, in Sec. [V] we end 
the work with some comments and conclusions. 



II. THE NEW PROCEDURE 

Let us start following the lines of ■ We consider the 
model described by a single real scalar field 0, with the 
non standard Lagrange density 



-X\X\ 



U(<f>) 



(1) 



2 



where 



and so we can write 



x = -d^d^ 

and U(4>) is the potential that identifies the theory. In 
this work we deal with bidimcnsional spacctimc, with 
metric ds 2 = dt 2 — dx 2 , using xo = t and a; 1 = — xi = x. 
Here we take h = c = 1 and we assume that the field and 
coordinates are all dimensionless. 

The equation of motion for this theory is 



->ra-l; 



d fi (\x\ n - 1 d»<t>) + u <l> = o 



(2) 



where \J $ = dU/d(f> and the energy-momentum tensor is 
given by 

T^ v = -g^C + Cxd^<f>d v (/), (3) 

where Cx = dC/dX. Here we are interested in static 
solutions, 4> = 4>{x), so we have that 



Too = P 
T n = 



1 



,/2n 



2n 

(2n-l) /2r 
2n 9 



U, 



U{4>). 



(4) 
(5) 



Moreover, the equation of motion ([2]) becomes 

(0' 2 "- 1 )' = U^ =<t>s ( x ). (6) 

where \^= ( p s indicates that we have to consider the field 
static, that is, <p = <Pi x )- This fact will be denoted from 
now on by \ s . This equation can be integrated once and 
we obtain 



In 



(7) 



where C is a constant that can be identified with the 
stress tensor T\\. Stability of the static solution imposes 
that C = 0, and this makes the energy density Too to get 
to the form 



P (x) = 4>'- 



p(x) 



2n 



2n- 1 



U(4>(x)). 



(8a) 



(8b) 



Let us now follow the procedure introduced in [18| to 
write the equation (|8aj) in another form, much more con- 
venient to study braneword models. The key step is to 
introduce a new function, W = W((f>), such that 



This fact leads us to 



(9) 



(10) 



0,n — 1 2n 



(11) 



Therefore, the energy density has the form 

p(x) = wf^ (12) 

We illustrate this procedure introducing two choices 
for the function W (</>). First, we choose the following 
n— dependent function 



( ^-2n+l;|0 2 



(13) 



this is a polynomial function of degree 4n — 1. The pres- 
ence of the hipcrgeometric function 2F1 introduces a new 
and general form to write W((f>); for example, we can 
write 



W(4>) = c^- l -<f 

W{4>) = 0-03 + |^-I^ 



(14a) 
(14b) 



for n = 1 and n = 2, respectively. 

In Ref. [17j . it was shown that the two choices given 
by Eqs. (fl"4"]) lead us with kinklike solution. Here we 
investigated the behavior of the model characterized by 
the general parameter n. 

Using the Eq. (fTTj) . we obtain the potential 



2n - 1 
2n 



1 



. 2n 



(15) 



For the first potential the equation (|10l) can be write as 
6' = 1 — 6 2 , whose the solution is 



4>{x) = tanh(x). 
The energy density is given by 

p(x) = scch 4n (x) 



(16) 



(17) 



The solution despite not having dependency of the pa- 
rameter, the energy density depends. The behavior 
asymptotic for different values for n is 



p(x) = Te 



+ 



(18) 



Is is clear that the solutions and energy densities asymp- 
totic behavior slow down with an increasing n. The en- 
ergy is 



E 



v^r(2n) 
T(2n+i) 



(19) 



For n = 1 and n = 2, the energy is E = 4/3 and E 
32/35, respectively. 

Now we introduce the second function, given by 



(20) 
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For example, we can write this function 



1 



W{4>) = |^T^^+-aresin(<£) 



(21a) 



W{4>) = fy/F^lQ - ^ 2 )+^arcsin(0)(21b) 

for n ~ 1 and n — 2, respectively. 

In Ref. [l7|, we investigated the case n = 2. Here we 
use the Eq. and we obtain the general potential 



u{4>) 



2n- 1 
2n 



1 



(22) 



The choice (f20|) leads us to 0' = ^/ 1 1 — 2 1 , which does 
not depend on n. It supports the compact solution 



for x < -f 
«/>(x) = <( sin(a:) for - f < x < f 
-1 for 5G > § 

with the respective energy density 

, for x < -§ 
cos 2n (x) , for -f < x < I 




(23) 



p(x) 



(24) 



, for x > ~ 



Here we note that the solution and energy density have 
a compact structure for all n. These kind of structures 
have also been studied in Refs [2l|, [22| . The energy of the 
solution is 



E 



^r(n + i) 

r(n+l) 



(25) 



For n = 1 and n = 2, the energy is _E = n/2 and -E 1 
37r/4, respectively. 



A. Linear Stability 

Let us now investigate linear stability. We introduce 
a small fluctuation rj(x,t) about the static solution <f>(x), 
that is, we write 



(x,t) = 4>{x) +r)(x,t) 



(26) 



where <f>(x) is solution of the static equation ([6]). With 
this, we obtain, at first order in 77, 



i/2n-2 



5^7?-2(n- 1)9^0 



d v (j)d v, q 



-U^T] = 0. 



Since <fi = <fr(x) is static solution, we can assume that 
r)(x,t) = T] a (x) cos(wt). Thus, we have 

- (2n - lM^-Y.]' + ^| s % = w 2 ^' 2 "- 2 7 ?s (27) 

We follow Ref. [l^, to rewrite the above equation as a 
Schroedinger-like equation. To do this, we introduce 



(z) = (2n - 1) V (n_1) Vs (V2n-l z) 



(28) 



which allows writing (|27)) as 

— u 22 (z) + u(z)u(z) = cj 2 u(z) (29) 

where 

nU^ n(n-l) 



V[Z) & n ~ 2 (2n- l)^™- 2 ' 
Then, using Eqs. (fT0|) and (|11| we obtain 



(30) 



2ti-3 
2n — 1 



v(z) = nW 4> 2n ^W^- 



2n- 1 



WL (31) 



For kinklike solutions given by Eq. ([13|) the potential v(z) 
is 



v(z) = 4(2n - l)n 2 - 2n(4n 2 - l)scch 2 (V2?i - 1 z) (32) 

This is the modified Poschl- Teller potential [H|. This 
potential supports the zero mode and other 2n — 1 bound 
states, with energies Ek = (2n — l)fc(4n — k), for k 
i). I ..... 2// I . All the others states of the model, with 
w > An 2 are not bounded. 

It is interesting to see that we could assoociate the 
parameter n with the number of bound states of the 
Schrocdingre-likc potential. As one knows, the number 
of bound states could in principe affect the rate of loss 
of the energy by radiation in dynamical processes, as for 
example, in a kink-antikink collision in comparation with 
the usual (j) 4 model (which is obtained with n = 1). See, 



e.g., Ref. [24 



For the compacton solutions given by Eq. (|2"0"|) . we get 

00 for z < —5^ 

v(z) =nA^-n+(n-l)sec 2 (Xz) for - ^ <z< ^ 
00 for z > 2j 

(33) 

where A = \/2n — 1. This is the Poschl- Teller potential 
[23I ]. The interesting feature of this potential is that it 
only supports bound states, and for k = 0,1,2,..., the 
corresponding eigenvalues are given by Ek = 4(2n — 
l)k(n + k). The radiation of energy for a collision be- 
tween compactons is then completely different from the 
case of kinks. 



III. TWINLIKE MODELS 

Let us now introduce a new family of twinlike models. 
Wc first recall that twinlike models are distinct mod- 
els having the same solution, and the very same energy 
density. The main objective here is then to introduce a 
family of models which is twin to the family of models 
investigated in the previous section. We follow the lines 
of Ref. @ and consider the theory 



C = -U(cf>)F(Y), 
where Y is defined as 

2"" 1 AIAI™- 1 

Y = 



n U{4>) 



(34) 
(35) 
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We note that if n = 1, we get back to the theory de- 
fined in also, for F(Y) = 1 + Y we obtain the model 
introduced in Eq. (JJ) above. 

This new model has the following equation of motion 

nd„ (lI^Fyd^ - U^F + YFyU$ = 0, (36) 

and the energy-momentum tensor is given by 

V = g^U(4>)F{Y) - nU^Fyd^d^. (37) 

where F Y = dF/dY. 

As before, here we are interested in static field config- 
urations; so, the equation of motion becomes 

-2n(u^F Y ^j + U 4> F -YF Y U^ = 0. (38) 

We suppose that Vi (i = 1, 2, n) is a set of static and 
uniform solutions of the equation of motion, meaning 
that U'(vi) has to vanish. Also, we use the energy density 
and take U(vi) — to make the energy itself vanish, for 
the static and uniform solutions. Recall that the same 
conditions work for the standard model. 

Since we are considering a new family of models, we 
guide ourselves with the null energy condition (NEC), 
that is, we take T^n^n u > 0, where n M is a null vector, 
obeying g^ u n tl n u = 0. This restriction leads to Fy > 0, 
for the general field configuration 4>{x, t) which is sup- 
posed to solve the equation of motion ([3"6")l . Moreover, 
for static solutions, the energy-momentum tensor gives 



Too - UF, 

Tn = -UF + 2nYF Y U. 
The Eq. ([38|) can be integrated once to give 

C 



2nYF Y - F 



U 



(39a) 
(39b) 



(40) 



Here we note that the solution <p(x) of this equation is 
the same solution <ft s (x) of the Eq. ([7]) , which appears for 
the previous model, with the position changed as x — > 
y/mx, with m — [c(2n — l)] 1 /". This means that we can 
write 

</>(x) = 4> s {^x), (44) 

and now the thickness of the solution is given by 

6 = S„/y/m. (45) 

Thus, the solution is thicker or thinner, depending on 
the value of c being lesser or greater than unit. We also 
note that c cannot be negative; and more, only stressless 
solutions have the specific form, given by Eq. (|44p . 

The energy density of the stressless solution ([4*3")) gets 
to the form 

p{x) = F{c)U{cl>{x)) (46) 
The energy is then E = F(c) dx U(<j>(x)), or better, 
F(c) 



E 



[c(2n - l)]V(3n) 
(2n - l)F(c) 
2n[c(2n- l)] 1 /^) 



E, 



dyU{Uv)) (47) 
(48) 



where E s is the energy for F = 1 + Y. 

The solutions with a non- vanishing Tu are different 
from the corresponding solutions of the previous model, 
because they do not have the form given by Eq. (|44p . 
We recall that for Tu = C, only the stressless solutions 
are stable [TU. Usually, the energy of the other possible 
solutions are divergent, and the solutions have oscillatory 
or divergent profiles. We find the same behavior in the 
standard model. 

Now we use again the formalism introduced in [l6l Kl\ 
to rewrite the energy density for the generalized model. 
Assuming that the Eq. ([9]) is valid, we have that 



Again, C is a constant identified with the stress tensor 
Tu. Furthermore, we have 



F v 



The potential is given by 



(49) 



Y 



1 



2n U(4>) 

The Eq. (|40|) can be written in the form 

C 



(41) 



(42) 



where G is a function with inverse G~ 1 (Y) = 2nYFy — F. 

For stressless solutions, that is, for C = 0, we have 
that 2nYFy=F and if we assume that G(0) = c, with c 
constant, real, we find that Y = c. With this result, we 
can rewrite Eq. (|4"2"]l in the form 



4>' 2n = 2ncU{(j>) 



(43) 



U{4>) 



F 2„-l 

2nY 



and we obtain the energy density in the form 



p(x) 



2nY 



-wi 



Now, using (|40|) with C = we have 
p(x) 



F Y ^W^-\ 



For m = 1, i. e., 



2n 



(50) 



(51) 



(52) 



(53) 



■5 



we have to impose 



Here we get 



Fy((2n-l)- 1 ) = l 



(54a) 



in order to make the Eqs. (|49l) and (|52j) identical to the 
Eqs. (fTUf and (fl2|). respectively. This also imposes that 



2n 



2n- 1 



(54b) 



The Eqs. (|54aj) and (|54b|) are the general restrictions on 
F(Y), to make the model twin of the previous model. 



A. Linear Stability 

Let us again investigate linear stability by introducing 
small fluctuations ri(x,t) in the static solution ^(a;), as 
done in Sec. Ill Al Using ([2"rJ|) in ([5rJ|) we obtain 



-2{n-l)F Y d^(j) 



-Fyd^rj + FyyAYd^cj) 



U^F-YFy)^ 
+U^{F-YF Y ) Y AY, (55) 



where 



Taking rj(t, x) — r) s (x) cos(wi) we get 

- [q{x) [2nF YY Y + (2n - 1)F Y ] v' s ]' = 



(56) 



U^(YF Y - F) 



//2n-l 



F YY Y 



u 



u 2 

p y2 u <i> 
-F YY Y — 



r) s + LU 2 F Y q(x)r] s 



(57) 



where q(x) = 0' 2n ~ 2 . 

In the case of a stressless solution we can use the 
Eq. (|40]) with C = to transform ((571) m the f° rm 



[g(x)^] ; + U H Yr, a = -^q{x)r,„. (58) 



where 



.4 



2 _ 2?iFy Y y + (2n - l)Fy 



(59) 



We note that A is constant for a stressless solution. Also, 
it is required that A is positive, if we want to ensure 
hyperbolicity of the differential equation. 

Again, we introduce the suggested exchange of vari- 
ables. 



u{z) = F$A*(l/ n - l Ti a (Az) 



(60) 



Uzz(z) + v 2 (z)u(z) = uj 2 u(z), 



(61) 



where 



v 2 (z) = nA 2 Y 



U, 



00 



l2n-2 



( n _l)y v 



i=<Ps{Az) 

(62) 



Now, using (j49|) and (|50| we obtain 

n(n-2) 



. ,4 2 F y ^ / W, 

«2(z) = — z : I n 



2n-l 



fef 2n-l 



(63) 



We note that if we impose the twin conditions (|54|). we 
obtain A 2 = 2n - 1 + 2n(2n - iy 1 F YY . With this, we 
obtain the relation V2(x) = v(x), if we choose 



F YY ((2n-l)- 1 ) = Q 



(64) 



Thus, we see that it is possible to find twinlike models 
starting from non standard theories. This is a new result, 
since before one usually started from a standard field 
theory, in order to construct the related twinlikc model. 



B. Examples 

Let us now specify the function F{Y) in order to il- 
lustrate how the formalism introduced above works. The 
first model we consider is 



F(Y) =a + bY\Y\ k 



k-1 



(65) 



Here we suppose that k > 1 and a, b are real number. 
To consider the stressless solution, we write 



2nkb - b 



i/fc 



(66) 



We note that 



F Y (c) = bk 



( a 



, , . 2nka 

and ^ (C) - 2^-T (67) 



\2nkb- b 
The conditions (|54)l give 
2nk - 1 

a ~ k(2n- 1) ' " k 
The function given by Eq. (|65p can be written as 



= (2n 1)k 1 (68) 



F(Y) 



2nk-\ (2n-l) fe " 1 
k(2n- 1) + k 



Y\Y 



fc-i 



(69) 



and so we get to the Lagrange density 



„ , N , .2^-^ X\X\ nk - 1 2nk-\ , 



(70) 
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To investigate linear stability we have to consider 

A 2 = 2nk - 1 (71) 

and so we note that the condition requires that k = 
1. However, in this case the twin theory is identically to 
the original model. 

We now introduce the second family of models, which 
obeys the strong condition. Let us consider the following 
function 



F(Y) =1+Y+G[Y 



1 



In- 1 



(72) 



For the three conditions (same solution, same energy den- 
sity and same stability behavior) to be valid, one imposes 
that 



G(0) = G'(0) = G"(0) = 0. 



(73) 



We can write a general function which obeys these three 
conditions; it has the form 



F ( y) = i + y + 5>(y-^-L_ 

i>2 v 



where all the fii are real parameters. 



(74) 



IV. BRANEWORLD MODELS 

Let us now investigate how the twinlike models studied 
above can be used to represent generalized braneworld 
models. Here we follow the lines of Ref. [l8|. In this 
context, we consider an action in five dimensions that 
describe gravity coupled to a scalar field in the form 



1 



S = / tPxy/g 1--R + £(</>, X) 



where we are using AttG = 1 and also 



x = \v M <t> 



M 



(75) 



(76) 



with M, N = 0,1, 2, 3, 4 running on the five-dimensional 
spacetime. The equation of motion which we obtain is 
given by 

G nm V n Vm4> + 2XCx4. - ^ = 0, (77) 

where 

G NM = C xx V M 4>V N <t> + g MN C X - (78) 

The energy-momentum tensor Tmn has the form 

Tmn = -QmnC + £xV 'm4>V ' n4> (79) 

The line element of the five-dimensional spacetime can be 
written as ds 2 = e^rj^dx 11 dx" — dy 2 , where A is used to 



describe the warp factor. We suppose that both A and 
cj> are static, such that they only depend on the extra 
dimension y, that is, A = A(y) and 4> = 4>(y). In this 
case, the equation of motion for the scalar field reduces 
to 

(2XC XX + Cx)4>" - {2XCx* ~ £*) = -^x4>'A! (80) 
Moreover, from the Einstein equations we get 



A" = -XC X 

A' 2 = l i {C-2XC x ) 



(81a) 
(81b) 



where X = —<fi' 2 /2 for static configuration, as before. 
To get to the first-order framework, we suppose that 



A' = -\w{<j>) 



(82) 



In this case, the Eqs. (|81a[) and (|81b[) lead us to, respec- 
tively, 



b'Cx = \w^, 



C-2XC X = 3^(0). 

In the case of theory (JT]) the equation of motion 
becomes 



(83) 
(84) 



(2n- 1 
We use Eq. ([531) to write 



/2n-2 ill , a j./2n— 1 a 1 



A' = U A 



2—w2 



(85) 



(86) 



and the potential 



*LjrWf* - V(0), (87) 



and the energy density 



bo = e 2A I 2T^W;- 1 - -W(0) 1 (88) 



1 1*7-2 I 

3 



Now we have to find the twin model. For this, let us 
consider a scalar field theory governed by the following 
Lagrange density 



C = -U(<l>)F(Y) + f(<j>) 



(89) 



where Y was defined in ([33]) and /(</>) is to be determined. 
We use the Eqs. (|53")l and to write, respectively 



(2F Y y- 



and 



b ' 2n = jr U ^> 



(90) 



(91) 



7 



where we have used f(<fi) = W 2 /3. 
We can so to rewrite (1911) in the form 



F = 2nYF 



Y 



(92) 



The Lagrange density of the twin brane model then 
has the following form 



C = -U(4>)F(Y) + ^W 2 (<j>) 



Moreover, the energy density is 



T, 



oo = e 



2.4 



771217-1 O 



(93) 



(94) 



which exactly reproduces the previous expression (|88[) if 
Fy(c) = 1 (95a) 

and, consequently, 

F(c) = 2n(2n- (95b) 

Thus, the two models have the same solution, with the 
very same energy density. And these are the two condi- 
tions required for the models to be twinlike models. 



and 



^fPayidahtf-dphaf,) = \w^d^4> (ioo) 



3 23 



-w^wh-^j^-w^ (ioi) 



The equation of motion for the scalar field leads to 



e 2A D4> - (2n - 1) 



W. 



2n-2 
2n-l It 



4n-3 



2^r 

4n— 3 

6 2 2 ™- 1 



(102) 



A. Brane Stability 

The investigation of the linear stability of the 
braneworld model can be done following Ref. [l8j]. The 
metric is perturbed in the form 

ds 2 = e 2A ^ ( Vlul + h^{y, x)) dx»dx v - dy 2 (96) 

and the scalar field in the form 

4> = 4>(y) + $(y,x) (97) 

For the starting model, given by Eq. ([1]), the first order 
contributions to the energy-momentum tensor are 
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The first order contributions to the Einstein equations 
are 
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For the general model (|93j) . after substituting the two 
twin conditions f|95[) . one is led to following set of equa- 
tions: i) the energy-momentum components: 
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(103b) 
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ii) the Einstein equations: 
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and 
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' 3 V 2n - 1 
iii) and the scalar field equation 
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We see that only for 
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(107) 



the set of equations is equivalent to that corresponding 
to the starting model. As we know, the study of stability 
is not a trivial task [25| ; however, we can assure here 
that, using the three conditions - same solution, same 
energy density, and the strong condition (|107p - the linear 
stability of the two models are the same. 

In the gravity sector, we can simplify the investigation 
of stability considering the transverse traceless compo- 
nents of metric fluctuations 

hfj, v = ^-(7r MQ 7r, y/3 + irppTTna) — -"K^a^j h al3 (108) 

where 7r M „ = rj^ — <9 M d„/D. Indeed, we can check that 



Eq. ([99]) reduces to the known equation 
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The next steps are known: wc introduce the z-coordinatc 
in order to make the metric conformally flat, with dz — 
e~ A ( y >dy and we write 



(110) 



In this case, the 4-dimensional components of h^ v obey 
the Klein-Gordon equation and the metric fluctuations of 
the brane solution lead to the Schrodinger-like equation 



where 



[-dl + U(z)} = p 2 H^ 



U(z) = \A'\z) + lA"{z). 



(Ill) 
(112) 



We note that the stability behavior in the gravity sector 
only depends on the warp factor A, so the two first condi- 
tions to make the two models twins - same solution and 
same energy density - are necessary for the two models 
to have the same stability behavior in the gravity sector. 

Therefore, we can write the following two impor- 
tant conclusions concerning stability of the two general 
models, described by Eq. ^ and by Eq. (|59")l . in the 
braneworld context: i) stability in the gravity sector is 
controlled by the warp factor, so the two first conditions 
for the models to be twins - explicitly, same solution and 
same energy density - lead the two twin models with the 
very same stability behavior; ii) stability in the scalar 
field sector are in general different, but the strong condi- 
tion - as given by Eq. (|107[) - makes the two models to 
have the very same stability behavior in the scalar field 
sector too. 

The above results show that the modifications pro- 
posed in the current work are robust and may be of direct 
interest to high energy physics. 



V. CONCLUSIONS 

In this work we introduced another route to construct 
twinlikc models, now starting from a generalized model, 
and generating another generalized model. We did this 
investigating several examples, showing that the proce- 
dure is generic and work for a diversity of models. 

To make the investigation stronger, in this work we also 
discussed the case of branes with warped geometry, in the 
scenario with a single extra dimension of infinite extent. 
Here we also investigated how the two first conditions for 
the models to be twins, namely, the same solution and the 
same energy density, and the extra condition, which we 
called strong condition, enter the game when one investi- 
gates stability. The result is that stability in the gravity 
sector is controlled by the warp factor, so it requires that 
the models are twins, that is, that they present the same 
solution, with the same energy density. In the scalar field 
sector, however, stability has also the same behavior if we 
includes the third condition, the strong condition (11071) . 

The procedure seems to be robust, working for several 
distinct models, valid both in the flat and curved space- 
time, in the last case for a braneworld model with a single 
extra dimension of infinite extent. 
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